Is the Zee model neutrino mass matrix ruled out? 
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Abstract 



A very economic model of generating small neutrino masses is the Zee model. 
£^ . This model has been studied extensively in the literature with most of the 



studies concentrated on the simplest version of the model where all diago- 



Oh, nal entries in the mass matrix are zero. SNO, and KamLAND data disfavor 



this simple version, but only when one also combines information from atmo- 



spheric and K2K data, can one rule out this model with high confidence level. 
We show that the simplest version of Zee model is ruled out at 3a level. The 
original Zee model, however, contains more than enough freedom to satisfy 
constraints from data. We propose a new form of mass matrix by naturalness 
consideration. This new form of mass matrix predicts that m U3 = 0, and 
tan 2 9 so i ar increases with \V e s\. For the best fit value of tan 2 9 so i ar , \V e %\ is 
sizeable but below the upper bound. 
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There are abundant data [1-6] from solar, atmospheric, laboratory and recent long base- 
line (K2K and KamLAND) experiments on neutrino mass and mixing. It is certain that 
some of the neutrinos have non-zero masses and also different neutrino spices mix with each 
other. In the minimal Standard Model (SM) in which there is just one Higgs doublet in the 
scalar sector and there are no right-handed neutrinos, neutrinos are massless. In order to 
have non-zero neutrino masses and mixing, one must go beyond the minimal SM. 

There are different possible ways to generate neutrino masses. A very economic way of 
generating neutrino masses is to introduce a charged scalar and an additional Higgs doublet 
into the minimal SM as proposed by Zee [7]. The Zee model provides a natural mechanism 
to generate small neutrino masses because they can only be induced at loop level, and also 
suggests special forms for the mass matrix. If one imposes a discrete symmetry such that only 
one of the Higgs doublets couples to the leptons as suggested by Wolfenstein [8] , one obtains 
a simple mass matrix with all diagonal entries zero. We will refer this simple version as the 
Zee- Wolfenstein model. This model has been studied extensively in the literature [7-12]. In 
this paper we further study the Zee model using the most recent experimental data. We 
show that the Zee- Wolfenstein model is ruled out at the 99.73% (3a) C.L.. However the 
original Zee model contains more than enough freedom to satisfy experimental constraints. 
We propose a new form of neutrino mass matrix resulting from naturalness condition. This 
model predicts that m vz = 0, and tan 2 9 so i ar increases with |V e3 |. For the best fit value of 
0.4 for tan 2 9 so i ar , \V e3 \ is sizeable but below the 3a upper bound. 

The Zee model contains, in addition to the gauge bosons and the minimal fermion con- 
tents, a singlet scalar h and two Higgs doublets 0i )2 transforming under the SM gauge group 
SU(3) C x SU(2) L x U(1) Y as (1,1,1) and (1,2,-1/2). With these particles it is not possible 
to have tree level neutrino masses from renormalizeable Lagrangian, but it is possible at one 
loop level. The relevant terms in the Lagrangian are [7], 

L = -T dR f*^i L e i3 - ^lf ab C^ bL e tJ h - WVpjPffrh + H.C., (1) 

where ip l aL = (v a L, &aLj and l a R are the left- and right-handed leptons with "a" the generation 
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index and the SU(2) L indices, e^- is the anti-symmetric symbol. C is the Dirac charge 
conjugation matrix. f^> ab are the Yukawa couplings responsible for charged lepton masses. 
f ab is an anti-symmetric matrix in generation indices a and b due to Fermi statistics. 

The mass matrix m for the charged leptons is given by, m = (viff + v 2 f 2 ) — v(sm Pff + 
cos P j "2). Here t> 7 =< 7 > are the vacuum expectation values (VEV), v = + v 2 = 174 
GeV and tan/5 = Vi/v 2 . The mass matrix m can be diagonalized to obtain the eigen-mass 
matrix m = Diag(m e , m M , m T ) by a bi-unitary transformation multiplying two unitary 
matrices Vl,r from left and right, m = V R mVl. 

The linear combination 0^, = cos/30]/ + sin/302 is "eaten" by W~. The physical 
combination which mixes with h is <p~ = sin/50]/ — cos/30^. We indicate the two mass 
eigenstates of masses Mi and M 2 for the charged scalars as hi = cos 9 z h — sin 9 Z (j) + and 
h 2 = sm9zh + cos9z4> + - Here sin6*z is proportional to M af3 characterizing the strength of 
the h — + mixing. 

The terms responsible for neutrino mass generation in the previous equation, in the mass 
eigenstates basis for the charged lepton and scalar fields, can be written as 

L = -E R mE L - E R { 1 m ^—f$)v L (sm9 z h\ - cos9 z h\) 

v tan p sin p 

-2ulfCE L (cos 9 z hi + sin 9 z h 2 ) + ... (2) 
where # = (f^ ab ) = V R f*Vl f = (f ab ) = VffVl E L , R = (e, ^ t) l , r , and v L = 

{yi,v 2 ,v z ) L . 

Exchange of charged scalars h 1>2 and charged leptons at one loop level, Majorana neutrino 
mass term L rn = (l/2)v'[M l/ Ci' L can be generated with 

M v = A[(fm 2 + m 2 f T ) - —Ajmjt + ffmf)}, (3) 

cos p 

where A = sm(29 z )log(Ml/Mf)/(16n 2 vtanf3) which is of order O(10" 5 ) if the sin(2^ z ) and 
tan (3 are both of order one. This is the general mass matrix in the Zee model [12]. The 
mixing matrix is the unitary matrix V which diagonalizes the mass matrix and is defined 
by, D = V T M U V, with D = diag(m Ul ,m U2 ,m U3 ). 
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The present experimental data on neutrino masses and mixing angles can be summarized 
as follow [13,14]. The 3a allowed ranges for the mass-squared differences are constrained to 
be: 1.6 x 10" 3 eV 2 < |Am 2 t J < 4.8 x 1(T 3 eV 2 , and 4.7 x 1(T 5 eV 2 < Am 2 solar < 1.7 x 10~ 4 
eV 2 , with the best fit values given by | Am 2 atm \ = 2.5x 10~ 3 eV 2 , and Am 2 solar = 7.0x 10~ 5 eV 2 . 
The mixing angles are in the ranges of 0.3 < sin 2 9 atm < 0.7 and 0.29 < tan 2 9 so i ar < 0.63. 
Also the CHOOZ experiment [4] gives an upper bound of 0.22 on the v e — v x (where v x can 
be either or u T or a linear combination) oscillation parameter for Am^ = Im^ 2 — \m Ul \ 2 > 
10~ 3 eV 2 . 

In the model discussed here the atmospheric neutrino and K2K data can be explained by 
oscillation between the muon and the tauon neutrinos, and the solar neutrino and KamLAND 
data explained by oscillation between the electron and muon (or a linear combination of 
muon and tauon neutrino) neutrinos. In this case the CHOOZ limit applies to the oscillation 
between the electron and tauon neutrinos *. 

Setting f$ in eq. (3) to zero, one obtains the famous Zee-Wolfenstein mass matrix, 



M„ 



a b 



a c 



y b C j 



(4) 



where a = Af^fa 2 — m 2 ), b = Af eT (m 2 — m 2 ) and c = Af^ T {m 2 — m 2 ). One can redefine 
the neutrino and charged lepton phases such that all a, b and c are real. 

Unfortunately the Zee-Wolfenstein model is now ruled out by experimental data. This 
can be seen from the following. 

The above mass matrix satisfies the "zero sum" condition m Ul +m„ 2 +m U3 = 0, therefore 
all the neutrino masses are determined in terms of the mass-squared differences [16]. We 
have [16] 



*There are additional evidences for oscillation between electron and muon neutrinos from LSND 
experiment [15]. If confirmed more neutrinos are needed to explain all the data. 
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(5) 



The other two masses are given by, m 2 2 = Am^ + m 2 Ui and m 2 3 = Am| 2 + m 2 U2 . The "zero 
sum" condition admits two types of mass hierarchy if the absolute value of r = Am^/Am^ is 
much smaller than one (experimentally |r| < 0.106 at 3a level), with one of them the normal 
one: m„ s > m V2 > m Ul and m Ul pa m„ 2 , and another the inverted one: \m V2 \ > \m Ul \ > \m V3 \ 
and m U2 pa —m Ul . One finds that \x\ = \m Ul /m U2 \ is determined to be very close to one. 

The mass matrix element M n = leads to, V 2 x m Vl + V 2 2 m V2 + V 2 z m vz = 0, which can be 
rewritten as 



Since \x\ is smaller but close to one, the above equation only allows negative x for small 
V 2 3 implying that only the inverted mass hierarchy is possible. One thus obtains a minimal 

v e2,min of V e2 close to (1 ~ V e3,max) / 2 ~ °- 47 > while data from SN0 and KamLAND prefers 
a smaller V 2 2 . Therefore SNO and KamLAND data disfavor the Zee-Wolfenstein model. 
This has been noticed in Ref. [11]. However, we would like to point out that although the 
Zee-Wolfenstein model can not produce the central values for the mixing and mass difference 
from solar and KamLAND data, the present data can not rule out the model at more than 
even 2a level. 

To have a more quantitative statement, we have carried out a detailed study and the 
results are shown in Figure 1. The dashed lines in Figure 1 are for tan 2 8 so i ar (sin 2 28 so i ar = 
4| V e i\ 2 \ K2I 2 ) with two values (0.22 and 0.15) of V e3 as a function of r. When |T4 3 | decreases, 
tan 2 ^ ; ar increases, tan 2 ^^ is about 0.53 for the 3a upper bound of \V e3 \, and becomes 
larger than the 3a allowed value of 0.63 when \V e3 \ decreases to be lower than 0.15. One 
therefore can take jV^I to be larger than 0.15 at 3a level. It is clear that the model is not 
possible to produce the best fit value of 0.4 for tan 2 9 so i ar - However at 2a, tan 2 9 so i ar can be 
as large as 0.54 [14]. Therefore it is not possible to rule out the model at more than 2a level 
from data on solar and KamLAND. 




x + (l + 2x)V 2 3 
1 — x 



(6) 
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FIG. 1. The dashed lines 51 and 52 are for tan Q solar &s functions of r. The two solutions for 
sin 2 a tm are indicated by solid lines Ala and A2a, and Alb and A2b, respectively. Here the indices 
"1" and "2" indicate the cases with \V e3 \ equals to 0.22 and 0.15, respectively. 

Data on sin 2 9 atm can provide further constraints on the model. The condition M 2 2 = 
m i^i + m 2^ 2 + m 3^3 — in the model can be used to determine sin 2 6 atm = V^ 3 . The 
mixing matrix V can be parameterized using three rotation angles, for example [1] V e2 = 
s i2Ci3, V e3 = s i3 and V^ 3 = S23 c i3- Here s^- = sinOy and = cos^-. Two of the angles, 
#12,13 can be determined in terms of V e3 and r from previous discussions. The condition 

M 22 = m Ul (s 12 + c 12 s 13 t 2 3) 2 + m U2 (c 12 - Si 2 si3^23) 2 + rn U3 cl 3 tl 3 = 0, (7) 

then determines 9 23 in terms of V e3 and r. Here £23 = s 23 /c 23 = tan^ 2 3- There are two 
solutions for tan 6 23 for given V e3 and y which we indicate by "a" and "b" . 

In Figure 1 the solid lines show sin 2 9 atm as a function of r for jV^I equals to its 3cr 
allowed upper value of 0.22 and the allowed lower value of 0.15. From the figure we see that 
sin 2 9 at m decreases for solution "a" , and sin 2 6 atrn increases for solution "b" as r increases 
from the 3a lower bound of -0.106 to the allowed upper bound of 0. All solutions for sin 2 9 atm 
are outside the 3a allowed range of 0.3 ~ 0.7. For |V e3 | smaller than 0.15, it is possible for 
sin 2 6 atm of solution "b" to become smaller than the 3a allowed upper bound. However \V e3 \ 
smaller than 0.15 will drive tan 2 9 so i ar to move out the 3a allowed range. Therefore the 
combined neutrino data on tan 2 Q solar and sin 2 9 atm rule out the Zee-Wolfenstein model at 
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more 3a level. 

The above discussions show clearly that the Zee-Wolfenstein neutrino mass matrix is in 
trouble. That does not, however, mean that the Zee model itself is in trouble. The mass 
matrix given in eq. (3) contains more than enough freedom to fit data. Here we encounter 
a common problem for physics beyond the SM that there are too many new parameters. 
Additional theoretical considerations have to be applied to narrow down the parameters. 

We find that an interesting neutrino mass matrix emerges if one requires that there 
should be no large hierarchies among the new couplings, that is all f tj and f 2 ah are of the 
same order of magnitude, respectively. This can be considered as a naturalness requirement. 
From eq. (3) one sees that all terms in the mass matrix are either proportional to m ; or mf. 
Since m T » m fte , the leading contributions to the neutrino mass matrix are proportional 
to f'^m 2 and f2 ,ab m T . To this order we have 

M u = -2A^—rft e m T , M 22 = -2-^-f^ft»m T: M 33 = 0, 
cos p cos p 

M l2 = —Mr It™ + F T ft e )m T , 

cos p 

M 13 = Arm T (m T - ^—ftl, M 23 = Aj^m T {m T - ^—ftl- (8) 

cos p cos p 

Without loss of generality, by appropriate choices of neutrino filed phases, the 11, 13, 23 
entries can be made real with just one physical phase S in the mass matrix. One can rewrite 
the above mass matrix as 



M v = a 



1 1 (ye i5 + x)/2 z ' 

(ye i5 + x)/2 xye iS xz 

y Z XZ j 



(9) 



with a = |M n |, x = \f» T \/\f e % y = \M 22 \/xa, z = \M 13 \/a. 

This is a highly constrained form of mass matrix. This matrix is rank two implying that 
one of the neutrinos has zero mass. The non-zero eigenvalues are given by 

2 

ml = —(^l + 2xycos5 + x 2 + y 2 ±^/{l + x 2 ){l + y 2 + Az 2 )) 2 . (10) 
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Since experimentally Am^ > 0, there are two types of eigen-mass hierarchies, a) m Ul = 0, 
\m U2 1 = \J Am 21 = m_, \m v . A \ = \J Am 2 2 — Am^ = m + ; and b) \m Vl \ = Am 2 2 | — Am^ = 
m_, \m U2 \ = Am 2 2 | = m + , m U3 = 0. The five parameters in the mass matrix are severely 
constrained from data on Am 21i32 , Ve2, V e % and V^. 

To have some idea about what parameter space may satisfy experimental constraints, 
let us discuss the situation with the phase 8 set to be zero for simplicity. For type a) of 
mass hierarchy since |r| = jAmf^/Am^l is much smaller than 1, one would have (1 + xy) 2 
to be almost equal to (1 + x 2 )(l + y 2 + Az 2 ). To satisfy this, x should be close to y and z 
to be much smaller than 1. Expanding the mixing matrix elements around x = y and small 
z, we find that (V e 2, V^, K-2) to be proportional to (z, xz, —(1 + x 2 )). Since z is much 
smaller than 1, one would obtain too small a V e2 in contradiction with solar and KamLAND 
data. This qualitative feature is not changed even if a non-zero 5 is introduced. There is no 
solution for the normal mass hierarchy of type a). 

For type b) of mass hierarchy, one has (V e3 V^ 3 , V r3 ) is proportional to (— 2xz, 2z,x — y). 
A small \r\ requires xy to be close to -1. Then small V e 3, and large jV^I and \V T 3\ require 
x to be small and 2xz to be of order one. We indeed find solutions for the mixing matrix 
satisfying experimental constraints. We also find that the size of V e3 anti-correlates with 
tan 2 9 so iar strongly, that is, when \V e ^\ decreases, tan 2 9 so i ar increases. If tan 2 9 so i ar is close 
to its best fit value of 0.4, \V e3 \ is close to, but below, the 3a upper bound of 0.22. In the 
following we present a sample solution with Am 2l 32 have their best fit values, 



m vi = 4.93 x 10" 2 eV, m U2 = -5.00 x 10" 2 eV, m„ 3 = 0. 

(11) 



^0.8244 -0.5312 -0.1953 ^ 



V 



0.2961 0.6989 -0.6511 



y 0.4823 0.4789 0.7335 j 



The tan 2 9 so i ar is 0.415 close to the best fit value. The value —0.1953 for V e3 is below, but 
close to the 3a allowed upper bound. 

In the above solution, the input parameters are: x = —0.3, y = 3.455, z = 1.667, 
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a = 1.94 x 1CT 2 eV. One can choose different signs for the parameters x, y and z. As long 
as the signs for x and y are simultaneously changed, the magnitudes of the eigen-masses 
and the mixing matrix elements are not changed. We will stick to the signs with x to be 
negative, y and z to be positive in our later discussions. 

One can also find solutions with smaller \V e3 \. For example, with x = —0.165, y = 6.531 
and z = 3.030, we obtain V e3 = —0.11, but tan 2 9 so i ar = 0.624 which is close to the 3a upper 
bound. 

We searched for other solutions. We find that it is also possible to have solutions with 
non-zero CP violating phase 5. For example with a = 1.92 x 10 _2 eV, x = —0.276, ye tS = 
3.467- i0.0573, and z = 1.571, we have 

m Ul = 4.93 x 10^ 2 e" il6 - 9 °eV, m U2 = -5.00 x 10-V n °eV, m V3 = 0. 

0.8147 -0.5048 — i0.2311 -0.1676 - i0.0024 ^ 



V 



0.3110 -i0.1995 0.7035 -0.6071 - i0.0087 

y 0.4166-i0.1619 0.4402 -^0.0561 0.7767 j 



(12) 



The value for tan 2 # so ; ar is about 0.464 which is within the la region. The value jV^I = 
0.168 is below the 3cr upper bound, but not far below. The CP violating Jarlskog parameter 
J = Im(ViiV 2 2V{ 2 V2i) is predicted to be —0.0165 which may be studied in future neutrino 
factories. We have kept masses in the form with phases to illustrate the existence of Majorana 
phases which can be rotated away by multiplying a phase matrix from the right on V 
obtained above. 

The neutrino masses obtained in the model are in the interesting ranges. The sum of 
the absolute neutrino masses, m sum = \m vi \ + \m V2 \ + \m U3 \, in this model is around 0.1 eV 
which is several times smaller than the recent bound of 0.69 eV from WMAP [17], but can be 
probed in the near future by the PLANK experiment where the sensitivity on m sum can be as 
low as 0.03 eV. Laboratory neutrino mass experiments can also test the model. A non-zero 
value a = \m ee \ can induce neutrinoless double beta decays. \m ee \ obtained here is about 
0.02 eV which is safely below the present bound [1,18] of 0.4 eV. However it can be probed 



by future experiments, such as GENIUS, MOON and CUORE, where sensitivity of about 
0.01 eV may be reached. The effective mass < m e >= \J\m Ul V e i\ 2 + \m U2 V e 2\ 2 + \m V3 V e3 \ 2 
measured by the end point spectrum of beta decay in our case is around ~ 0.05 eV which is 
unfortunately a factor of 2 smaller than the sensitivity of future KATRIN experiment. 

The off-diagonal entries of the couplings f ab and /.f'" 6 can induce flavor changing inter- 
actions. One should make sure that constraints on related parameters will not rule out the 
regions of the parameters to produce the mass matrix discussed above. It is not possible to 
completely determine the couplings using just information from neutrino masses and mix- 
ing. We therefore take a simple situation with f2 ,TT = for illustration. In this case for the 
example given in eq. (12): f$' Te / cos/? = -0.33 x 10 -2 , ff ™/ cos/? = -(1.14-i0.02) x 10 -2 , 
Af er = 0.93 x lO-^GeV- 1 ), and Af» T = -0.26 x 10- n (GeV -1 ). It is interesting to note 
that the solution obtained here is consistent with the naturalness requirement that f2 ,Te to 
be the same order of magnitude as f$ ' Tfl , and f er to be the same order of magnitude as f^ T . 
If one chooses a smaller x one would obtain bigger hierarchy for the parameters, and 
f er . The qualitative features will not change when other values for the parameters are used. 

Exchange of the neutral Higgs boson fa (with mass M ) can induce at tree level k — > IjlJk 
decays. For the values of f$ ,Tfl and /vf' 1 " 6 obtained in the example of eq.(12) we have 

B(t -> wfi,nee) w 3.5 x lO -9 ^, 0.80 x lO -13 ^; B(t -> /ry) « 0.76 x lO -8 ^; 
B(t -> efjLji,eee) « 2.9 x lO" 10 ^, 0.67 x lO" 14 ^; B(r -> e 7 ) « 0.63 x lO" 9 ^. 

In the above B T = (100(GeV)/M tanf3) A B SM (r -> i/ T ^P„) with B sm (t -> v T ^w^ « 17%. 

There are experimental constraints on the above decays with the 90% C.L. bounds given 
by [1]: B(t -> ^ft, fxee) = 1.9 x lO" 6 , 1.7 x 10~ 6 , B(r -> e/^f, eee) = 1.8 x lO" 6 , 2.9 x lO" 6 , 
S(r — > [/,"/, cy) = 1.1 x 10~ 6 , 2.7 x 10 -6 . For tan/? of order one, all the branching ratios 
predicted above are safely below the experimental values if the mass M is of order 100 GeV. 

Non-zero / u can also induce radiative charged lepton decays by exchanging charged 
scalars. If the parameter A is not too much smaller than a natural value of A = 10 -5 
(GeV -1 ), their contributions for the rare decays mentioned will be much smaller. The rare 

10 



decays of the types discussed in the above will not provide significant constraints. 

From the above discussions we see that the new form of mass matrix proposed is con- 
sistent with present experimental data. It also predicts m„ 3 = and a sizeable | V^3 1 . In 
particular, if the error on tan 2 9 so i ar is reduced and the present best fit value holds, | V^3 1 will 
be close to the 3a allowed upper bound. The model can be tested in the future. 
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